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Abstract 

In this paper the existence of a class of self-similar solutions of the 
Einstein- Vlasov system is proved. The initial data for these solutions are 
not smooth, with their particle density being supported in a submani- 
fold of codimension one. They can be thought of as intermediate between 
smooth solutions of the Einstein- Vlasov system and dust. The motivation 
for studying them is to obtain insights into possible violation of weak cos- 
mic censorship by solutions of the Einstein- Vlasov system. By assuming a 
suitable form of the unknowns it is shown that the existence question can 
be reduced to that of the existence of a certain type of solution of a four- 
dimensional system of ordinary differential equations depending on two 
parameters. This solution starts at a particular point Po and converges to 
a stationary solution Pi as the independent variable tends to infinity. The 
existence proof is based on a shooting argument and involves relating the 
dynamics of solutions of the four-dimensional system to that of solutions 
of certain two- and three-dimensional systems obtained from it by limiting 
processes. 



1 INTRODUCTION 



It is well known that solutions of the Einstein equations coupled with suitable 
models of matter can yield singularities in finite time. The unknowns in these 
equations are the spacetime metric and some matter fields. The exact nature 
of the latter depends on the physical situation being considered. The usual 
terminology in general relativity is that there is said to be a singularity if the 
metric fails to be causally geodesically complete, i.e. if there are timelike or 
null geodesies which in at least one direction are inextendible and of finite affine 
length. The singularity is said to be in the future or the past according to the 
incomplete direction of the geodesies. It is expected on the basis of physical 
intuition, and known to be true in some simple cases, that the geodesic in- 
completeness is associated with the energy density or some curvature invariants 
blowing up. For background on this subject see textbooks such as [T3], [H] and 
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[2"5] . One of the best known types of singularities in general relativity are those 
which occur inside black holes. When a black hole is formed by the collapse of 
matter it is known that under suitable circumstances an event horizon is formed 
which ensures that the singularity can have no influence on distant observers. 

Mathematical relativity is the study of the properties of solutions of the Ein- 
stein equations coupled to various matter equations. One of the main questions 
in the field is the cosmic censorship hypothesis. There are two versions of this 
conjecture called weak and strong cosmic censorship, both of which were pro- 
posed by Roger Penrose. It should be noted that, contrary to what the names 
might suggest, the strong version does not imply the weak one. The results 
proved in what follows are motivated by weak cosmic censorship and strong 
cosmic censorship will not be discussed further here. Weak cosmic censorship is 
a statement which concerns isolated systems in general relativity. Mathemat- 
ically this means considering solutions of the Einstein equations which evolve 
from asymptotically flat initial data. Initial data for the Einstein equations con- 
sist of a Riemannian metric h a b, a symmetric tensor k a b and some matter fields 
which for the moment will be denoted generically by Fq, all defined on a three- 
dimensional manifold S. Solving the Cauchy problem for the Einstein-matter 
equations means embedding the manifold S into a four-dimensional manifold 
M on which are defined a Lorentzian metric g a p and matter fields F such that 
h a b and k a b are the pull-backs to S of the induced metric and second funda- 
mental form of the image of the embedding of S while Fo is the pullback of 
the matter fields. The metric g a p and the matter fields F are required to sat- 
isfy the Einstein-matter equations. A comprehensive treatment of the Cauchy 
problem for the Einstein equations can be found in [26]. Initial data on R 3 are 
called asymptotically flat if the metric h a b tends to the flat metric at infinity 
in a suitable sense while k a b and F tend to zero. Physically this corresponds 
to concentrating attention on a particular physical system while ignoring the 
influence of the rest of the universe. 

A solution of the Einstein-matter equations evolving from initial data is said 
to be a development of that data if each inextendible causal curve intersects 
the initial hypersurface precisely once. When this property holds the initial 
hypersurface is said to be a Cauchy hypersurface for that solution. In general, 
a solution is called globally hyperbolic if it admits a Cauchy hypersurface. For 
prescribed data there is a development which is maximal in the sense that any 
other development can be embedded into it. It is unique up to a diffeomorphism 
which preserves the initial hypersurface. 

In a spacetime evolving from asymptotically flat data it is often possible 
to define future null infinity I + as a set of ideal endpoints of complete future- 
directed null geodesies. We can say that any singularity occurring does not 
influence events near infinity if there is no inextendible causal curve to the fu- 
ture of the initial hypersurface which is incomplete in the past while intersecting 
a future-complete null geodesic. The first of these properties means intuitively 
that this curve represents a signal which comes out of a singularity while the 
second property means that it reaches a region which can communicate with 
infinity. If a curve of this type does exist it is said that a globally naked singular- 
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ity exists. The past of null infinity, J~(X + ), is the set of points for which there 
is a future-directed causal curve starting there and going to null infinity. The 
complement of J~(l + ) is called the black hole region. Its boundary is called 
the event horizon and is a null hypersurface in M. 

There is a notion of completeness of null infinity. A precise definition will 
not be given here but roughly speaking it corresponds to the situation where 
there are timelike curves contained in J~(X + ) which exist for a infinite time 
towards the future. Physically this means that there are observers which can 
remain outside the black hole for an unlimited amount of time. If the maximal 
globally hyperbolic development of asymptotically flat initial data always has a 
complete null infinity then this ensures the absence of globally naked singulari- 
ties. For any inextendible causal curve to the future of the initial surface which 
goes to null infinity must intersect the initial hypersurface. Hence it cannot be 
incomplete in the past. The completeness of X + ensures that the solution is 
large enough to represent the whole future of a system evolving from the initial 
data under consideration. The intuitive content of the weak cosmic censorship 
hypothesis is that in the time evolution corresponding to initial data for the Ein- 
stein equations coupled to reasonable (non-pathological) matter the existence 
of a singularity implies that of an event horizon which covers the singularity 
and hides it from distant observers. Often this is weakened to the requirement 
that a horizon exists in the case of generic initial data. Up to now this intuitive 
picture has only been developed into a precise mathematical formulation under 
special circumstances. In general finding the correct formulation is part of the 
problem to be solved. 

Due to the mathematical complexity of the Einstein equations many of the 
studies related to singularity formation for these equations have been carried 
out for spherically symmetric solutions. In spherical symmetry the Einstein 
vacuum equations are non-dynamical due to Birkhoff's theorem, which says 
that any spherically symmetric vacuum solution is locally isometric to the 
Schwarzschild solution and, in particular, static. Thus it is essential to in- 
clude matter of some kind. A matter model which has proved very useful for 
this task is the scalar field. This is a real-valued function <fi which satisfies the 
wave equation V Q V Q = 0. In this case the Einstein equations take the form 
Ra/3 — 87tVq,</)V^0, where R a p is the Ricci curvature of g a p. The spherically 
symmetric Einstein-scalar field equations were studied in great detail in a series 
of papers by Demetrios Christodoulou. This culminated in [7] and [8]. In [7 
it was shown that in this system naked singularities can evolve from regular 
asymptotically flat initial data. This represents a problem for the weak cosmic 
censorship hypothesis but the conjecture can be saved by a genericity assump- 
tion since it was shown in [5] that generic initial data do not lead to naked 
singularities. 

For the spherically symmetric Einstein-scalar field equations it is known 
from the work of Christodoulou [B] that small asymptotically flat initial data 
lead to a solution which is geodesically complete and hence free of singularities. 
(In fact this small data result has recently been extended to the case without 
symmetry [15].) On the other hand there are certain large initial data for which 
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it is known that a black hole is formed. The threshold between these two types 
of behaviour was studied in influential work by Choptuik [3] and many other 
papers since. This area of research is known as critical collapse and is surveyed 
in [TJ] . It is entirely numerical and heuristic and unfortunately mathematically 
rigorous results are not yet available. 

The scalar field provides a simple and well-behaved matter model. At the 
same time no such field has been experimentally observed and the matter fields 
of importance for applications to astrophysics are of other kinds. One astro- 
physically relevant matter field which has good mathematical properties is col- 
lisionless matter described by the Vlasov equation. The necessary definitions 
are given in the next section. For the moment let it just be noted that the un- 
known in the Vlasov equation is a non-negative real- valued function f(t, x a , v b ) 
depending on local coordinates (t,x a ) on M and velocity variables v b . Ana- 
logues of a number of the results proved for the scalar field have been proved 
for the Einstein- Vlasov system. For small initial data the solutions are geodesi- 
cally complete [22] ■ There are certain large initial data for which a black hole is 
formed [T] . The threshold between these two types of behaviour has been inves- 
tigated numerically in [23] and [18]. A closely related matter model which has 
been very popular in theoretical general relativity is dust, a fluid with vanishing 
pressure. It is equivalent to consider distributional solutions of the Vlasov equa- 
tion of the form f(t,x a ,v b ) — p(t,x a )S(v b — u b (t,x a )) where the 5 is a Dirac 
distribution. From many points of view dust is relatively simple to analyse. 
Unfortunately it has a strong tendency to form singularities where the energy 
density blows up, even in the absence of gravity. For this reason it must be 
regarded as pathological and of limited appropriateness for the investigation of 
cosmic censorship. A detailed mathematical study of formation of singularities 
in the Einstein equations coupled to dust was given in [5]. In spherical sym- 
metry dust particles move as spherical shells. It can easily happen that shells 
including a strictly positive total mass come together at one radius and this 
causes the density to blow up. This effect is known as shell-crossing. 

The motivation for this paper is the wish to understand cosmic censorship 
better for spherically symmetric solutions of the Einstein- Vlasov system. Is it 
true that in asymptotically flat spherically symmetric solutions of the Einstein- 
Vlasov system there are no naked singularities for generic data so that colli- 
sionless matter is as well-behaved as the scalar field? Could it even be that the 
Vlasov equation is better-behaved and that there are no naked singularities at 
all? No answers to these questions, positive or negative, are available although 
considerable effort has been invested into obtaining a positive answer. In what 
follows we try to obtain new insights by approaching a negative result through 
an interpolation between dust and smooth solutions of the Vlasov equation and 
looking for self-similar solutions. There are some results on related equations 
which give some hints. In the case of the Vlasov-Poisson system, the non- 
relativistic analogue of the Einstein- Vlasov system, global existence for general 
data, not necessarily symmetric, was proved by Pfaffelmoser [20] and Lions and 
Perthame [16]. The relativistic Vlasov-Poisson system, which is in some sense 
intermediate between the Vlasov-Poisson and Einstein- Vlasov systems, (but not 
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in all ways) has been shown to have solutions which develop singularities in fi- 
nite time. Rather precise information is available about the nature of these 
singularities [13]. 

As a side remark, we mention a paper [27] where it was suggested that 
naked singularities are formed in solutions of the Einstein- Vlasov system. The 
solutions concerned were axially symmetric but not spherically symmetric. The 
work is purely numerical but trying to understand what it means for the an- 
alytical problem leads to the conclusion that the solutions computed in [37] 
were dust solutions rather than smooth solutions of the Einstein- Vlasov system. 
This is discussed in [23] . There are also reasons for doubting that the numerical 
results really show the formation of a naked singularity [2J5] . 

A class of distributional solutions of the Einstein- Vlasov system intermedi- 
ate between smooth solutions and dust is given by the Einstein clusters [TT] . 
These are spherically symmetric and static, i.e. there exists a timelike Killing 
vector field which is orthogonal to spacelike hypersurfaces. It is supposed that 
the support of / consists of v a such that the geodesies with these initial data 
are tangent to the spheres of constant distance from the centre of symmetry 
on these spacelike hypersurfaces. This means that the radial velocity and its 
time derivative in the geodesic equation are zero. These are in general two inde- 
pendent conditions on the data at a given time. A wider class, the generalized 
Einstein clusters [5], [3J, is obtained as follows. In the case of the Einstein clus- 
ters taking the union of the spheres at a fixed distance from the centre defines 
a foliation of the spacetime by timelike hypersurfaces and the condition on the 
support means that the four-velocity of a particle with the given initial data 
is everywhere tangent to these timelike hypersurfaces. The generalized Ein- 
stein clusters are obtained by dropping the condition of staticity and replacing 
the family of timelike hypersurfaces invariant under the timelike Killing vector 
field by another foliation by timelike hypersurfaces which intersect any Cauchy 
surface in spheres and whose equation of motion follows from the Vlasov equa- 
tion. Once again the four- velocity of a particle in the support of / is tangent 
to these hypersurfaces at all times. An analytical formulation of this definition 
will be given in the next section. It should be noted that the generalized Ein- 
stein clusters exhibit shell-crossing singularities and thus can still be thought 
of as pathological. We are interested in them as an intermediate step towards 
better-behaved matter models. 

There are two major differences between the generalized Einstein clusters 
and the solutions studied in this paper. In the case of Einstein clusters the 
value of the angular momentum of the particles F is uniquely determined by 
the distance r to the centre of symmetry. By contrast, in the solutions studied 
in this paper the angular momentum takes a continuous range of values for each 
value of r. The second difference is that in the case of Einstein clusters at each 
spacetime point the component of the velocity vector v a of a particle in the 
direction of the vector d r takes on only one value. In the case of the solutions 
obtained in this paper the component along the direction d r of the velocity 
vector takes on two different values at most spacetime points. This only fails 
at some exceptional values of r at a given time. This difference in the structure 



5 



of the generalized Einstein clusters and the solutions considered in this paper is 
what gives some plausibility to the idea that the solutions described here could 
be a big step towards better-behaved matter models. From the physical point 
of view, in the case of the generalized Einstein clusters the material particle 
with the smallest value of r would not experience any gravitational field, and 
therefore could not approach the centre r = unless its angular momentum 
vanished. In the solutions studied in this paper, since two radial velocities are 
allowed at each spacetime point, the material particle with the smallest value 
of r changes in time. This allows the occurrence of a collective collapse of the 
whole distribution of particles towards the origin with some of them coming 
closer and closer to the center as the value of some suitable time coordinate t 
increases. 

Self-similar solutions of the massless Einstein- Vlasov system have also been 
considered in the paper |17j . There are several differences between the approach 
in [T7] and the one considered in this paper. The first one is the choice of the 
rescaling group under which the solutions are invariant. The massless Einstein- 
Vlasov system is invariant under a two-dimensional group of rescalings. The 
choice of a particular one-dimensional rescaling group has been made in this 
paper by imposing that the distribution function / for the particles remains 
always of order one (see Sections [U . This condition is natural, because the 
function / is invariant along characteristic curves. On the contrary, the choice of 
one-dimensional rescaling group for the solutions in 17 imposes that / becomes 
unbounded near the singularity for the particles within the self-similar region, 
something that can be achieved assuming that the distribution of matter is 
singular near the light-cone. The second difference between the solutions in [17] 
and those in this paper is that the solutions in [17] can be thought of as self- 
similar perturbations of the flat Minkowski space. As a matter of fact they have 
been computed by means of a perturbative iteration procedure that takes flat 
space as a starting point and where the terms in the resulting series have been 
computed numerically. By contrast, the solutions of this paper are obtained by 
means of a shooting procedure in which a parameter that measures the amount 
of energy in the self-similar region is of order one. The approach in this paper 
uses purely analytical methods and docs not rely on numerical computations. 
On the other hand, in order to simplify the arguments, we have restricted the 
analysis in this paper to the study of dust-like solutions, an assumption that 
was not made in [17] . 

The plan of the paper is as follows. We will first reduce the problem of 
finding self-similar solutions of the Einstein- Vlasov system to an ODE problem 
that can be transformed into a four-dimensional system using suitable changes 
of variables. Using these transformations it will be seen that the construction 
of the desired self-similar solutions reduces to finding a particular orbit in the 
corresponding four-dimensional space connecting a certain point with a steady 
state that has a three-dimensional stable manifold. The existence of such an 
orbit will be shown by adjusting a parameter that measures the density of par- 
ticles in a particular perturbative limit. The precise limit under consideration, 
which has the goal of making the problem feasible using analytical methods, cor- 
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responds to assuming that the radius of the region empty of particles, measured 
in the natural self-similar variables, is small. 



2 THE EINSTEIN-VLASOV SYSTEM IN 
SCHWARZSCHILD COORDINATES. 



We do not use exactly the classical Schwarzschild coordinates, but a slight mod- 
ification of them that normalizes the time to be the proper time at the center 
r = 0. The metric is given by (cf. |21|): 



ds 2 



Mr) dt 2 + e 2Mt,r) dr 2 + ^ ^2 + ^2 ^2) 



(2.1) 



If we restrict our attention to spherically symmetric solutions it is convenient 
to use the quantities (cf. [2"T]): 



ii X ' v |2 

r = \x\ , w = , F — \x A v\ 



to parametrize the velocity variables. In particular F is constant along charac- 
teristics. Writing the particle density as 

f = f(r,w,F,t) 

the Einstein- Vlasov system for spherically symmetric solutions in these coordi- 
nates becomes: 



fi—X 



F 



where: 



E 



1 + W 2 + -T 



and the functions A, /i that characterize the gravitational field satisfy: 

e~ 2X (2r\ r - 1) + 1 = 87rr 2 p, 
er 2X (2rfi r + 1) - 1 = 8nr 2 p 

with boundary conditions: 

M(0)=0 , A(0) = 0, 
A(oo) = 0. 

On the other hand p and p are given by: 



/oo " />oo 
/ EfdF 
-oo UO 

a f°° r r°° w 2 



dw. 



(2.2) 

(2.3) 

(2.4) 
(2.5) 

(2.6) 
(2.7) 

(2.8) 
(2.9) 
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With these basic equations in hand it is possible to give some details con- 
cerning generalized Einstein clusters, as promised in the introduction. These 
are not required to understand the main results of the paper but help to put 
those results into a wider context. A distributional solution of the Vlasov equa- 
tion whose support is a smooth submanifold £ has the property that S is a 
union of characteristics of the equation. A simple example is that of dust where 
the support is the graph of a function u a (t,x) of the form W(t,r)^-, When 
expressed in terms of polar coordinates this becomes the graph of a function 
W(t, r) augmented by the condition F = 0. Here the function W solves the 
equations 

dt E 

^- = -(X(t, R)W + e^^-^'^n'tt, R)E) (2.11) 
dt 

where E = Vl + W 2 . 

Now consider the generalized Einstein clusters. They are only defined under 
the condition of spherical symmetry. They can be thought of as defining a 
matter model which can be used in the spherically symmetric Einstein-matter 
equations. Here they will be described in terms of Schwarzschild coordinates. 
The basic unknown is a function R(t,r) which satisfies i?(0,r) = r. It is the 
area radius at time t of the shell which had area radius r at time 0. As input 
we require a function F(r) which is the angular momentum of the particles on 
the shell which was at radius r at time zero and N(r) which is the density of 
particles per shell evaluated on the shell which had area radius r at t = 0. For 
some purposes it is more convenient to use R as a radial coordinate instead of r 
and this is what was done in the original papers [9] and [2] . For a given shell at 
a given time the angular momentum and radial velocity of the particles are fixed 
and so the intersection of the support of the solution with the fibre of the mass 
shell over the point with coordinates (i, r) has codimension two. The following 
equations should be satisfied: 

d 4 = e^, (2.12) 
dt E v ' 

^-(XW + ^B-^^) (2.13) 



where E = J 1 + W 2 + and the functions A and /i are to be evaluated at the 
point (t, R). These are the full characteristic equations for the Vlasov equation. 
The difference in the coupled system comes from the fact that the expressions 
for the components of the energy-momentum tensor are different in the two 
cases. In the case of Einstein clusters the characteristics of interest have W = 
and =0. It follows immediately that = 0. In that case the angular 

momentum is related to the geometry by the relation F — r M 



\ — r\i' ' 

The equations which have been written up to now describe particles of unit 
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mass. We are interested in the construction of solutions of (|2.2[) - (|2.9p supported 
in a region where (w 2 + t) takes large values near the formation of the singu- 
larity. This suggests replacing (|2.3I) by: 



£ = ^ + _ (2 . 14) 

The system (|2~2|) . (|2~4l) - (|2~9| . (|2.14|) is invariant under the rescaling: 

r -^9r , t-*6t for t < , u> -> -4= to , F ^ OF (2.15) 

for any > 0. It is then natural to look for solutions of (|2.2p . (l2~4l) - (| 2. 14 1) 
invariant under the rescaling (|2.15j) . They will be the self-similar solutions in 
which we will be interested in this paper. 

The system obtained when (|2.3p is replaced by (|2.14p can be interpreted as 
describing particles of zero rest mass. The rationale for this assumption is that 
near the singularity the derived solution will satisfy w 2 + ^ >> 1, and therefore 
it could be expected that it is possible to treat the whole Einstein- Vlasov system 
with massive particles as a perturbation of the massless problem. 

In what follows we will consider solutions of (|2.2p . (I2.4I) - (|2.14I) where / is 
not a bounded function, but a measure concentrated on some hypersurfaces that 
will be described in detail later. As was mentioned in the introduction there 
is a class of distributional solutions of the Einstein- Vlasov system which are 
equivalent to what is usually known in the literature as dust. From this point of 
view the solutions considered in this paper are intermediate between dust and 
smooth solutions and hence will be called dust-like solutions. Note, however, 
that in contrast to dust they do have some velocity dispersion. The dimension 
of the support of / in the tangent space at a given spacetime point is zero for 
dust, one for generalized Einstein clusters, two for the solutions in this paper 
and three for smooth solutions. For the solutions here it will be possible to 
describe the distribution of velocities for the particles at a given point using a 
function depending on one coordinate, while a general distribution of velocities 
compatible with the assumption of spherical symmetry would depend on two 
coordinates. 



3 SELF-SIMILAR SOLUTIONS 



In this section we formulate the system of equations satisfied by the solutions 
of ([272]) . (f274|) - ([279l) . (f2~T4|) that are invariant under the transformation (f27l5|) . 
We will call these self-similar solutions in what follows. It is convenient, as a 
first step, in order to transform (|2.2p . (|2.4p - (|2.14p to a more convenient form to 
define a new variable: 

» = (3-D 

Vf 
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We will assume in the rest of the paper that / = for (r,v,F,t) — (r,v,0,t) 
in order to avoid singularities in (|3.ip . Moreover, we can even assume a more 
stringent condition on /, namely / = for < F < Sq for some So > 0. 
Concerning the support in the r coordinate, the solutions constructed in this 
paper will vanish for r < yo (— t) for some yo > 0. 

Making the change of variables (r, w,F,t) —> (r,v, F, t) and denoting the new 
distribution function by / with a slight abuse of notation we can transform the 
system ^T^, $Z ty - $rfty into: 



d t f + e»- x T d r f 
E 



X t v + e^HrE - e^ x —^ )d v f = 
r A E 



E 



4/ E 



fFdF 
fFdF 



dv, 
dv. 



(3.2) 

(3.3) 
(3.4) 

(3.5) 



Notice that the change of variables (|3.1|) eliminates the dependence on the vari- 
able F for the characteristic curves associated to the Vlasov equation (cf. (|3.2I) V 
Moreover, the functions p and p and therefore the functions A, /i characteriz- 
ing the gravitational fields depend on / only through the reduced distribution 
function: 



C(r,«,t) 



fFdF. 



(3.6) 



In particular, it is possible to write a closed problem for the reduced distribu- 
tion function that can be obtained multiplying p. 21) by F and integrating with 
respect to this variable: 







E M 


E = 




+ 1 




^2 J 


OO 


P = 


E(dv 7 

-oo 






—Qdv 

-oo E 


P = 



^- x p r E - e»- x 



1 



r^E 



d v ( = 0, 



(3.7) 

(3.8) 
(3.9) 
(3.10) 

The system ()3.7[) - (|3.10[) complemented with (|2.4I) . (|2.5p is a closed system of 
equations. 

We will now study the class of self-similar solutions of the system (|2.4p , (|2.5p . 
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(|3.2[) - (|3.5p . These are the functions having the functional dependence: 

f(r,v,F,t) = G{y,V,$) , n(r,t) = U(y) , X(r,t)=A(y), (3.11) 
r „ , , F 



y 



v = {-t)v , $ 



(3.12) 



(-*) v ' (-*) 

The solutions of (|2.4I) . (|2.5[) , (I3.2l) - (|3.5p with this functional dependence satisfy: 



2/G y - VG V + ^G* + e 



U-A 



V 

E 



G u 



yAyV + e u ' A U y E - e u ~ A 



y 3 E 



where: 



and 



E= JV 2 + — 



(2yA y - 1) + 1 = 8ny 2 p, 
(2yU y + 1) - 1 = 8ny 2 p 



-2A 



with boundary conditions: 



U = , A = at y = 0. 



Here: 



7T 

^ = ~2 

y 

7T 



G$d$ 



00 y2 p r°° 



00 -B Uo 



dV. 



(3.13) 
(3.14) 

(3.15) 
(3.16) 

(3.17) 

(3.18) 
(3.19) 



The function G which is a solution of (|3.13[) - (|3.19[) is constant along the 
characteristic curves of p,13|) which are given by: 



dy 
da 



U-A 



V 



v 2 + ± 



y + e 



u-A 



Vy 



dV 
da 

da 



-V- yA y V + 



v 

U - A u„ 



y 



WV 2 y 2 + l-e 



U-A 



(3.20) 



(3.21) 



(3.22) 



In these equations a is just a parameter that is used to parametrize the charac- 
teristic curves. Its precise definition will be given later in some specific cases. 
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The equations (I3.20l) - (|3.22l) can be integrated explicitly for any pair of func- 
tions U = U (y) , A = A(y). Indeed, the first two equations can be rewritten 
as: 

An, £>TJ 

(3.23) 
(3.24) 

uu uy 

where: 



dy _ 


A dH 


da 


dV 


dV _ 


_ e _ A dH 


da 


dy 



H = — y/V 2 y 2 + l + yVe A . (3.25) 
V 

The trajectories in the (y, y)-plane associated to the solutions of (|3.20[) . (|3.21[) 
are contained in the level sets: 

H = h. (3.26) 

We will also need the self-similar formulation of the integrated form of the 
equation ()3.7|) . In this case the function £ in (|3.6|) has the functional dependence: 



((r,v,t) = (-t) 2 e(y,V). 

Notice that: 

&(y,V) = I G$rf$. (3.27) 

The function 9 satisfies: 



o 



y&y - VBv - 29 + e u ~ A Z- Q v 

yk y V + e u - A U y E - e u ~ A -^r ) 9 y 

y 3 EJ 

= (3.28) 

and: 

/oo 
EedV, (3.29) 
- OO 

7T f°° V 2 

P=~2 —j-QdV. (3.30) 
The characteristic curves associated to p. 281) are p.20p . Q3.2ip and: 

de 



29. (3.31) 
da 
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4 SINGULAR SELF-SIMILAR SOLUTIONS: 
GENERAL PROPERTIES. 



The main goal of this paper is to construct a family of distributional solutions of 
(|3.13l) - (|3.19p for which G — G (y, V, <&) is a measure supported on some surfaces 
in the three-dimensional space with coordinates [y, V, $) . In this section we 
will describe in a heuristic manner the argument yielding the construction of 
such solutions. The arguments will be made rigorous in the rest of the paper. 
The key idea behind the argument is that the problem can be transformed into 
a system of ordinary differential equations for the particular class of solutions 
described in this section. 

Taking into account that the singularities of the distribution G might be 
expected to be propagated by characteristics it is natural to look for solutions 
of (13~T31 - (13~T31) of the form: 

G (y, V,$) = A (y, V, $) S (H (y, V) - h) (4.1) 

satisfying (|3 . 13[) in the sense of distributions. Let us assume that A, H have the 
differentiability properties required for all the following formal computations. 
Plugging flU} into (|3T3)l we obtain: 

(a (y, V) A y + b (y, V) A v + S (H - h) 

+ A(a (y, V) H y + b (y, V) H v ) 5' (H - h) 
= 



where: 



a( y ,V) = y + e u -^ = ^, (4.2) 
b (y,V) = -V-(yA y V + e u - A U y E-e u - A ^)=-e A ^. (4.3) 

V y EJ 



dy 



Notice that a (y, V) H y + b (y, V) H v = 0. Then: 

(a (y, V)A y +b (y, V) A v + $A $ ) S(H-h) = 0. 

This equation is satisfied if: 

a(y,V)A y + b{y,V)A v + <f>A* = (4.4) 

on the surface {H = h} x K+. Let us assume that the curve {H = h} can be 
parametrized, at least locally, using a parameter a satisfying: 

y = y (a) , V = V (a) , 

^P-=a(y(a),V(a)) , ™M = b (y (a) , V (<?)) . (4.5) 
da da 
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Then the function A can be written on the surface {H = h} x K + as a function 
of the variables (a, <£>) . We can write: 

A (y (a) , V (a) , $) = A (a, $) for (y, V, a) e {H = h] x M+ (4.6) 

and using (|4.5p we can rewrite (j4.4j) as: 

A CT + = 0. (4.7) 

Since the curves {H — h} can be determined in terms of alone it is con- 
venient to compute this distribution explicitly. If G has the form (|4.1[) the 
distribution defined in (|3.27|) is given by: 

@(y,V)=/3S(H-h) (4.8) 

where: 

poo 

/3= 

Since A is given by (|4.6I) it follows that: 

/•OO 

P = /3(a)= l(a,$)$d$ for (y, V) € {F = h} . (4.9) 
Jo 

We can compute /? (cr) along the curve {H = h} . To this end we multiply (|4.7I) 
by $ and integrate in the $ variable in the interval [0, oo) . Then: 

Pa - 2/3. 

The function /3 then takes the form: 

(3 {a) = /3 e 2ff (4.10) 

for some /3q > 0. 

In the rest of the paper we prove that there exist functions A (a, $) as in 
(glp and curves {# = h} with A, [/ solving ([535]) - ([5TT]) and p, p as in ([5TT5]) . 
(|3.19p such that (I4.1[) solves p.!3p in the sense of distributions. 



5 SINGULAR SELF- SIMILAR SOLUTIONS: 
DESCRIBING THEIR SUPPORT. 

In this section we describe in a precise manner the form of the curved surface 
containing the support of the distribution G for the self-similar solutions con- 
structed in this paper. Such a surface is contained in the surface S = 7 x M. + , 
where 7 C {(y, V) : y > , V £ R} is an unbounded curve, at a strictly positive 
distance from the line {y = 0} = {(y, V) : y — , Vel} with a discontinuity 
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in its curvature at the point (yo,Vo) £ 7 placed at the minimum distance from 
the line {y — 0} . In order to avoid such irregular curves it is more convenient 
to assume that the curve 7 is the union of two analytic curves 71 and 72 that 
can be parametrized in the form: 

~/i = {(y, V):y <y< 00, V = Vi(y)} , i = l,2 (5.1) 

where the functions Vi (y) are analytic and satisfy: 



lim Vi (y) = lim V 2 (y) = V = - 



Vi (y) < V 2 (y) for y Q < y < 00 



(5.2) 
(5.3) 



for some yo £ (0, 1) . Since the curves 7, are contained in the curve {if = h} it 
follows that the functions Vi (y) are the two roots of the equation: 



y 

assuming that such roots exist. Then 
1 



VvV + 1 + yVe A = h 



(5.4) 



Vi (y) = 



V2(y) 



2U 



y 2 e 2A ) 



ye A h - a / (ye A hy - (e 2U - y 2 e 2A ) ( — - h 2 



.21! 



(5.5) 



( e 2C/ _ y 2 e 2A) 



-ye A h + a / (ye A hf - (e 2U - y 2 e 2A ) 



-:,2U 



h 2 



(5.6) 



Notice that for such solutions the support of G in (|4.1j) is contained in the 
half-plane {y > y } ■ Therefore, p(y)=p (y) = for y < y . Then (|3~T5|) - ([3~T7| 
imply U (y) = A (y) = for y < y . 

Under suitable regularity assumptions for the curves 7^ near the point [yo, Vq) 
that will be made precise below the functions U and A are continuous at the 
point y — y^. In such a case (|5.4p implies: 



yo 



yo 



(5.7) 



We will prove later that it is possible to construct the desired curves 7, , i = 1, 2, 
defined by means of (|5.1[) with the property that the following limits exist: 



Vi (y) - V a 



lim 



Ki 



Ki £ 



1,2, K x <K 2 . 



(5.8) 
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Moreover, the quotients of the functions A and U by y/y — yo also tend to finite 
limits. Let: 



hm -#L = lS 
]im -^L = 9 2 e 

y^vt Vv - yo 



(5.9) 
(5.10) 



We parametrize the curve 7 = {H — h} as in the previous section using a 
parameter a. We will denote a parameter of this kind on the curves 71, 72 by 
ci, er 2 respectively. Due to (|4~2"j) , (|4~5]) . (1531) it follows that: 



dcji 



dy a(y,Vi(y)) y + e u - A , vdv)y 



, i = 1,2. 



(5.11) 



We will normalize the parameters oi = iji (y) by means of the condition: 

<7i(y )=0 , i-l,2. (5.12) 

Finally we remark that in order to obtain the functions U and A we need 
to prescribe the distribution 8 defined by (j3~2"7j) . Using P~g]) . (|4~TU)) it then 
follows that: 

BnYf ^ lP 2cr l(») 



Iff (vMv))] 

^X{y>y„}e 2 ^ y) 

Iff (v,v 2 (y))\ 



S(V-V 2 (y)) 



(5.13) 



where X{y>y } is the characteristic function of the half-plane {y > yo} . Using 
([339]) . (pHO]) it follows that: 



p(y) 



7T/3oX{y>yo} 



2<ri(y) 



o2cr 2 (t/) 



+ ■ 



If (l/.Vad/)) 

7r/?0X{y>y o } 



Iff (y, ^1(2/)) 



■\/(V 2 (y)) 2 y 2 + i 



■\/(Vi(y)) 2 y 2 + l 



(5.14) 



o2o- 2 (y) 



(Vi(y)Y 



Iff (»,Vi(y))| y/fy^ytfy* 



(V 2 (y)Y 



Iff y/(y 2 (y)) 2 y> + i 



(5.15) 



and the functions t/ and A can then be obtained using the equations (13 . 1 5[) . 
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Due to the dust-like character of the solutions considered in this paper, they 
exhibit a singular behaviour for p and p at the radius y — yo- This singularity 
is due to the fact that at this point the radial velocity of the particles, in self- 
similar variables, vanishes. However, since the motion of the trajectories after 
they reach the singularity continues in a smooth way, and since p and p are 
integrable near this radius, this singularity can be expected to disappear if the 
dust-like assumption is relaxed and some thickness is given to the support of 
the distribution function in the phase space. 

The main result of this paper is the following: 



Theorem 1 There exists £o > small such that, for any tjq G (0,£o) there exist 
a value of f3$ > and two curves 71 , 72 that can be parametrized as in (5.1]) 
with the functions Vi (y) , V% (y) as in (5.5]) . \5.6\) satisfying (5. 6 2\) . (5.3]) . \5.8\) . 
the functions U, A satisfying (3.15]) . (3.16]) and (5.9^) . (5.10]) with p, p as in 
\5.14\) , (5.15]) and ax, a 2 solving (5.11)) . (5. IS]) . 



Using Theorem[T]it is possible to obtain distributional solutions of the prob- 
lem (|3.13|) - (|3.19[) . In order to make the definition of the distribution G in 
(|4.1I) precise we use (|4.6I) . (|4.7[) . Let us prescribe a smooth function Aq ($) in 
$ G (0, 00) . Taking into account (|4.7I) we can then define: 

A(er,$) = A (e~ CT $) . 

Using the structure of the curves 71 , 72 it would then follow that the distribution 
G in (|4.ip would be given by: 



^0 ( e -ifa)q>) X{y >yo} 
A (e-«(v)$) X{y>yo} 



g (y, v, *) = v ) dH - J;™r* s ( v ~ Vi (y)) 



\av{y,V 2 (y)) 
We then have the following result: 



ohT_ TA S(V-V 2 (y)). (5.16) 



Theorem 2 Suppose that the function Aq (•) G Cq (0, oo) satisfies 

/ ^ ($)$d$ = (3 . (5.17) 
Jo 

Let us define a Radon measure G € M (K + x K x M + ) by means of (5.16)) with 
the functions V\ (•) , V 2 (•) , <J\ (•) , o~2 (•) as in Theorem^ Then the functions 
p, p defined 113.18]) . (3.19]) belong to the spaces Lf oc (0, oo) for 1 < p < 2. 
The functions A, U defined by means of (3.15]) - (3. 17\ ) belong to (0, oo) for 
1 < p < 2. T/ie measure G satisfies 113. 13]) in the sense of distributions. 
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Remark 3 The space Cq (0, oo) is the space of compactly supported continu- 
ously differ •entiahle functions and the space M (IR + x K x R + ) is the space of 
Radon measures on ]R + x M. x M+. It is not necessary to require Aq (•) to be 
compactly supported. Actually this condition could be replaced by assumptions 
of fast enough decay near the origin and infinity. 

Remark 4 It is worth noticing that the functions p, p associated to the distri- 
bution G have an integrable singularity as y — » y$ . 

In the rest of this section we will prove Theorem[5J Theorem [T] will be proved 
in the remaining sections of the paper using a shooting argument and refined 
asymptotics of the solutions for yo small. The following auxiliary result will be 
used in the proof of Theorem [2] and it will be proved in Section [B] We remark 
that Theorem [2] will not be used in either the proof of Theorem [T] or that of 
Proposition [5] below. 



Proposition 5 The curves 71, 72 whose existence has been proved in Theorem 
[I] satisfy the following conditions: 



lim $(v,Vi(y)) 



tim m^m 



= u 



for some constants L\ < L 2 . 

Proof of Theorem [H Using ([3^7]) . (|5~16|) and (fSTT]) we obtain: 



(5.18) 



9 (y, V) 



f3 e 



2<ri (y) 



X{j/>J/o} 



s(y-vi(y)) 



OV 



X{y>yo} 



Iff (y,v 2 (y))\ 



6(V-V 2 (y)) . 



(5.19) 



We can then compute p, p using (|3.29[) , (|3.30p 

nPoX{y> yo } 



p{y) = 



^(v)Jl + yi (Vi (y)f e 2 ^Jl + yi (V 2 (y)) 



P(y) = 



^PoX{y>y } 



Iff (y,Vi(y))\ 



e a«n(») ( Vl (y)Y 



Iff (2/^(2/)) | 



(5.20) 



dH 
L I dV 



(y,Vi(y))\ ^1 + yHV^y)) 2 



e 2°2(y) (y 2 {y)f 



'ff (y,V 2 (y))\ Jl + y* (V 2 (y)) 



(5.21) 
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Using (j5~^j) - ([5~TUj) . ([57X5]) . we obtain: 



\p(y)\ + \P(y)\ < 



CX{y>y } 

\[y - 2/0 



(5.22) 



whence the estimate p, p € Lf oc (0, oo) , 1 < p < 2, in the theorem follows. On 
the other hand ([XT5|) - ([XT7|) imply: 



A 



U 



1 1 A &7T 

log 1 

2 V y 



» [(8< 2 p(Q + l)e 2A ^ -1] 

n 2l 



(5.23) 
(5.24) 



Due to Theorem Q] the functions A, [7 are bounded for any finite value y > 0. 
On the other hand, (fB33"l) . (|5T24f imply A, [/ e W^'J (0, oo) , 1 < p < 2. 

In order to conclude the proof of Theorem [5] it only remains to prove that G 
solves (|3.13l) in the sense of distributions. This is equivalent to showing that: 



yAyV + e u ~ A U y E ~ e L 



y 3 E 



= o 



for any ip = tp (y, V, $) € Cg 



oo 



99 I GdydVd® 
' v- 

(5.25) 

~) . Using (15.16[) we can rewrite (|5.25l) 



as: 



i=1 Jy Q JQ 



oo poo 



{y^p) + (Vp) v - ($cp), 



E 



yA y V + e u - A U y E - e u - A -^r ) p 



1 

y J E 



A ( e - CT4 ^)$) 



(5.26) 



■d$dy 



= 



and making the change of variables e ai ^)$ — >• $ we can transform J into: 
2 



2/o ->o 



yAyV + e v - A U v E - e 



M-A 



1 



y 3 .E 



9 



($)e^ 



(y,Vi(y),*e< T i(!')) I 9V 



^(V,Vi(|/)) 



dyd& 
(5.27) 
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Notice that we can write: 



F = — ( W ) y + (Vlf) V ~ (<M $ - {f~ K \v 



yA y V + e u - A U y E - e u ~ A 



V 3 EJ v 



-yip y + V(p v - ($<^) $ - e 



u-A 



Vy 



and, using Leibniz's rule: 

F = -yA y ip - ye A (e~ A <p) y + Vip v - -f-U, 



A/. A „, !• . ;u - - tj e U-A V V 



M-A 



V 



--Lp + e 



U-A 



V 3 y 2 



V 1 + V 2 V 2 (1 + V 2 y 2 ) 



-w — e ; (e y>) 

v/TT¥v v ;y 



+ h/A y + CT v e 



(7-A 



e 



(7-A 



e u - A U v 



v/i + yy-' 



- (l + v^v) 5 
i 



(7-A 



After some cancellations: 



' "¥>)„- *e A ( e_i V)*-P 



+ ( y A y y e A + y e A + ^i^yr^Tv 



1 



(5.28) 



Then (15.271) can be rewritten as 

2 />oo />oo 



v v 1 \%v (y> v i(y))\ 



= o. 



Due to Proposition [5] as well as the fact that the curves 71 , 72 are globally 
dchncd it follows that: 



dH 
W 



(y,Vi(y)) 



4)1-1 W {y ' vdy)) 



Then: 



^=E(- 1 ) i " 1 / / r 



00 poo 



^(2/,V-(y)) 



(5.29) 
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Using (PT25|) and (j5~2gf we obtain: 
F (y, V, $) 



ye A + e U K£ 



y / V 2 y 2 + l 
1 



y + e 



u-A Vy 



(e"V) 



y + e 



fj-A 



1 



,£/-A 



Equations (jO|l . ([Q|) and (pTTTj) give: 



1 



y/V 2 y 2 + l 



yA y V + V+ ^v^TW 

V y 



,a-A 



(5.30) 



Therefore 



^(y,Vi(y)) 



i(y) 



da 



1 ( „-A, 



It then follows, using the chain rule that: 



V(y,Vi(j(),*e '<{«')) 



(e' f(v) e- A !"V (y, *5 (l/) , = -f (y, Vt (y) , $e^) ^ „ , ( ^ . 



Formula (|5.29p then becomes: 



J = £ (-I) 1 " 1 ^ A (*) A (e^^e- A ^V (y, V5 (y) , $e^)) d$dy 



or, equivalcntly: 



„- i ** 
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and using ([5^]) . (j5lU|) . (|5T2j) : 

2 />oo 

J = V (-I) 1 " 1 / 1 (*) V? {vo, Vi (y ) , $) d*. 

Due to the fact that ip (vq, V, (uq) , $) = ip (yo, Vq, $) for i = 1,2 we have J = 
and (|5.26[) follows. This concludes the proof of the theorem. ■ 

We now remark that it is possible to derive some detailed information about 
the behaviour of the curves 71, 72 as y —¥ 00. 



Theorem 6 Suppose that the curves 71 , 72 are as in Theorem [7J Then, the 
following asymptotic formulas hold: 



U = log \jL\ + log \Jl-y^ + o (1) as y 
A -> log (Vs) as y — > 00, 



00, 



7i = _ 2W3(l -^) (1 + 0(1)) „ 
(l-4^)y 

V2 = — _ — - — — (1 + 0(1)) as y —> 00 



V3yo y \ y 

for a suitable constant C\ G ffi. 

Notice that the asymptotic behaviour of the solutions in Theorem [5] shows 
that the support of these solutions approaches the line {V = 0} away from the 
self-similar region (i.e. for y — > 00). This is the one of the main differences 
between the solutions described in this paper and the ones in |17j . 

It is relevant to notice that the spacetime described by the solutions in 
Theorem^ exhibits curvature singularities and not just coordinate singularities. 
To this end we use Kretschmann scalar (cf. [25]): 

R a ^ s R aM s = 4K 2 + + l2r- 2 V a V b rV a V b r 

where K is the gaussian curvature of the quotient of the spacetime by the 
symmetry group and to is the Hawking mass that can be computed by means 
of: 

m = - (1 - d a rd a r) . 

Combining (I2.1j) . (I3.11[) , p,12[) we obtain the following self-similar form for the 
Hawking mass: 

to = - 1 — e v F*) 1 
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and therefore, it follows from Theorem [5] that: 



r r 
m ~ - for - — - sufficiently large. 
3 (-t) ' 5 

On the other hand, the last term in the Kretschmann scalar can be written as 
(cf. [10], Appendix A): 

24r- 2 (^(k- V h rV c r) + 27rrtrr) +96tt 2 ( T ab - ^gab) ( T ab tlT - ab 



The last term turns out to be positive for any matter model satisfying the dom- 
inant energy condition, which includes in particular the case of Vlasov matter. 
Therefore R a ^ 1& R a p 1 s > and so the curvature becomes singular as r — > 
for a fixed large value of . 

We remark that the solutions which have been derived do not provide an 
example of violation of the cosmic censorship hypothesis for Vlasov matter, be- 
cause the spacetimes concerned are not asymptotically flat asr-> oo. Moreover, 
it turns out that the region contained inside the light cone reaching the singular 
point at r = 0, t = 0~ in the spacetime described by Theorem[6]is dependent on 
the data on the whole region with < r < oo. This implies that a gluing of this 
spacetime with another one causally disconnected from the singular point is 
not possible, because this would require doing some gluing along regions where 
r = oo. In order to check these statements it is convenient to rewrite the metric 
(|2~Tj) in double null coordinates. Notice that ([2"7T]) . (|3~TT|) and (j3~T2]) yield the 
following self-similar structure for the metric: 

ds 2 = -e 2U (r^))dt 2 + e 2A (<^))dr 2 + r 2 (d9 2 + sin 2 Odip 2 ) . 

The double null coordinates are then just the constants of integration associated 
to the pair of differential equations: 

-e U ( T=% )dt + e A ^)dr = 0, 
e U ^)dt + e A ( T^))dr = 0. 

The solutions of these equations can be written in terms of two integration 
constants u and v that will define the double null coordinates. The particular 
choice of coordinates has been made in order to obtain u and v taking values in 
compact sets: 

rv e A(£)-U(£) 
arctanh (u) = log (-t) + - - ^ Aii) _ u{0 

arctanh (v) = log (-t) - - _ j^p^ d£ 

In the region close to the centre (i.e. y « I) the structure of the metric 
is similar to Minkowski. On the other hand, Theorem [6] yields the following 
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asymptotics for r >> (— t) 



arctanh (v) 



log(-t) + 



log 



( 



r 



) 



arctanh (u) 



log(-i) + 



log 



(-*) 



r 



) 



The light cone approaching the singular point is described in these coordinates 
by the line u = 1. Notice along such a line, for v of order one we would have 
r = oo, whence the assertion above follows. 

For these reasons a spacetime behaving asymptotically as Minkowski can- 
not be obtained gluing the self-similar solution obtained in this paper with a 
spacetime causally disconnected from the singular point. This kind of gluing 
might be possible for non self-similar solutions of the Einstein equations behav- 
ing asymptotically near the singular point like those described in this paper. 
However, such an analysis is beyond the scope of this paper. 



The strategy used to prove Theorem Q] is the following. We first transform 
the or iginal problem (I3T51) . (I3T6)) . (Oil . (IQl) . (1531) . (15^1) . (l5^1)- (l5T2l) . (I5T41) . 
(|5.15[) into a family of four-dimensional autonomous systems depending on the 
parameter /3q by means of a change of variables. It will be shown that proving 
Theorem [T] is equivalent to finding an orbit for this system connecting two 
specific points P , Pi of the four-dimensional phase space. The point Pi is 
a unstable saddle point with an associated three-dimensional stable manifold 
A4 = M. (/?o) that can be described in detail in the limit yo — > 0. A shooting 
argument will show that for a suitable choice of the parameter /3o the manifold 
Ai ((3q) contains the point Pq. In the rest of this section we give the details of 
this argument. 

6.1 Reduction of the problem to an autonomous system. 

Instead of the set of variables (y, U, A, Vi, <Ji) it is more convenient to use the 
set of variables (s, u, A, Q, Qi) where: 



6 PROOF OF THEOREM QJ 
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Then, the evolution equations (|3.15p . (|3.16[) , (|5.4p . (|5. become: 



dQi 
ds 



= <Jl-y* , i= 1,2 , Ci < C2, 



-2A 



(2A S - 1) + 1 



Qi [C1 2 + 1] 



1,2, 



01 [CI + 1] 



(6.2) 
(6.3) 



Cie" + y eVCi +1 C 2 e" + y e A y/@ + 1 



(6.4) 



-2A 



(2u s +3) - 1= - 



Qi (C1) 2 



>i(C 2 ) 2 



Cie" + y e A VCl+1 C 2 e" + y n e A VCl+1 



where 



16tt 2 /3o 

2/0 



The initial conditions p,17[) . (|5 . 1 2[) imply: 

u = , A = , Qi = l , i = 1,2 at s = 0. 



(6.5) 
(6.6) 

(6.7) 



Notice that the system (|6.3[) - (16.5[) with Q as in (|6.2p is a four-dimensional au- 
tonomous system of equations for the unknown functions (Qi, Q2, A, u) . Notice 
however that the system seems to becomes singular if the variables (Qi, Q2, A, u) 
approach the values in (|6.7p due to the vanishing of the denominators in 



(|6.5p . To treat these singularities we rewrite the terms 
Notice that (|6.2p implies: 



Jft>e A ^/(Ci) 2 + l + Cie u 



1 



(1 - y2 e 2(A-«)) 



-yoVl-% 2 e A - 2u T^ 



i = 1,2, 



^(e- 2 " (1 - y 2 ) - 1) (1 - y 2 e 2 ( A -)) + % 2 (1 - % 2 ) ^(^0. 



(6.8) 



(6.9) 



Then: 



yoe A V (0) + 1 + Cie" = Te"^ , i = 1, 2 
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and the system of equations (|6.3p - (|6.5p becomes: 
dQi Q1C1 



ds 
dQ 2 
ds 

- 2A {2A S - 1) + 1 



Z ' 
Q2C2 



6e~ 



-2A 



{2u s + 3) - 1 



f [C 2 + i]+f [C 2 2 + i] 



? 2 d 2 , QiCi 



z 



z 



(6.10) 
(6.11) 
(6.12) 

(6.13) 



We now eliminate the variables A, u in (|6.3p - (|6.5|) and replace them by the 
functions Z and G where Z is as in (I6.9[) and G is defined by means of: 



G 



-2A 



Then f|6. 12[) becomes: 

G s = 1 - G - 
On the other hand 



8e~ u 



'Ql 



-2u 



2 

implies: 

Z 2 + l 

[(i-vD + vh 



[d 2 + i]+$[C 2 2 + i] 



1)G 



[G + y 2 (l-G)] 



whence: 



u = — loe 
2 S 



(Z 2 + 1) G 



JG + % 2 (1-G)] i 
Differentiating this formula we obtain: 

ZZ S vl G fl 



s (Z 2 + l) 2 [G + y 2 (l-G)]G- 
Eliminating u s from this formula using (|6.13[) . (|6. 15|) we obtain: 

'3 1 



where: 



zz, = 



2G 



-A (Z 



1 



4G 



Q?(Ci) 2 , Ql (C2) 2 



2 [G + y 2 (l-G)]G 
Using (|6.15[) it then follows, after some computations, that: 
4GZ [G + yl (1 - G)] A = 2 (1 - G) j/ 2 Z 
+ 0e- [-y 2 [Ql [d 2 + 1] + Q\ [C2 2 + 1]] 
+ [Q 2 (Ci) 2 + Q2(C2) 2 1 [G + ylil-G)] 



Z 



(6.14) 



(6.15) 



(6.16) 



(6.17) 



(6.18) 
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The last bracket in (|6 . 18[) can be rewritten as: 

-i/o [Ql [C? + 1] + Ql [C 2 2 + 1]] + [0? (CO 2 + Ql (C2) 2 ] [g + yl (1 - G)] 
= Ql [Cl - vl (d 2 + 1)] + Ql [C 2 2 - vl (Ci + 1)] 
+ (i- % 2 ) [q? (Cif + Q 2 . (C2) 2 ] (G-i). 

Using (I6.8[) we obtain: 



(6.19) 



[C 2 -2/o 2 (C 2 + l)] 



(l-y 2 )Z 2 ^(l-^'e 



2\ 2 c A-2u 



± 



(1 -y 2 e 2 ( A -«)) 2 (l-y 2 e 2 ( A -")) 2 



% 2 



2^ e 2(A-2„) 



- 1 



i = 1,2. 



(1 - % 2 e2(A-)) 
Plugging (|6~20|) into (|6~T9|) it then follows that: 

-vl [Ql [Ci + 1] + Ql [C2 2 + 1]] + \q\ (Ci) 2 + Q2 (C 2 ) 2 1 [G + y 2 (1 - G)] 



(6.20) 



(1 -y 2 e 2 ( A -«))- 



{Ql 



2yo (1 - 



2 _A— 2tt 



+ 2/0 



(l^oT 



2(A-2u) 



(l-y2 e 2(A-«)y 

(Qi + Q 2 ) 



(Q 2 - Q 2 ) Z 



(1 - y 2 e 2 ( A -«)) 

+ (1-2/0) [Qi (Ci) 2 + q!(C 2 ) 2 ] (G-l) 

and using (|6.18l) we arrive at: 

A= {l-G)yl 

2G [G + yl (1 - G)\ 

0e~ u f (1 ~ Vl) Z 



4G [G + yl (1 - G)] [ (1 _ % 2 e 2(A-„)) 



2 (Qi + Ql) 



2yo(l 



2^2 e A-2u 



(l-y 2 e 2 ( A -«))' 



(Q? - Q 2 ) + ^ 



(6.21) 



where: 



2(A-2u) 



(Q? + Q 



_ (1 - ^e 2 ^-")) 
+ (l-y 2 ) \Ql (Ci) 2 + Ql (C 2 ) 2 1 {G-l) 



(6.22) 



In order to obtain analytic solutions it is convenient to introduce the change of 
variables: 

ds = 2GZdx , X = at s = 0. (6.23) 
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Then the system (f6"TT0|) . (|6~TT|) . (1(05)) . (|6~TT|) becomes: 
dQ x 



dQ 2 



2GQiCi, 
-2GQ 2 C 2 , 



= 2G 



Z(l-G)- 



dG 

dZ 
dx 

with the initial conditions: 



6e~ 



[Ql [Cl + ^+Ql [Cl + l]] 



^ = (3G - 1 - 2GA) (Z 2 + 1) 



Qi = Q2 = 1 , G = 1 , Z = , at x = 0. 
We can further simplify $ in (|6.22[) using 



2(A-2u) 



_ (1 - J/2 e 2(A-«)) 

(l-?/ 2 )(G-l) 
(l- % 2 e 2 ( A -")) 2 



1 



{Qi + Q 



(6.24) 
(6.25) 
(6.26) 
(6.27) 

(6.28) 
(6.29) 



(l-y 2 ) e 2 ^+Z 2 } (Ql + Ql) 



2yjl - y*Ze A - 2u (Qi - Q 



In order to identify the behaviour of $ as Z — > we write the terms in brackets 
on the right-hand side of (|6.29j) as: 



(1 



,2(A-2u) 



(1 - y 2 e 2 ( A -«)) 
1 

" (l-2/ 2 e 2 ( A -«)) 2 



Then becomes: 

vl (1 - vl) (Qi + Q 



1 



(l-% 2 ) 2 (e 2(A - 2u) 



1+2(1 



vl) vl ( ( 



2(A-u) _ 1 



$ = 



(1 - t/ 2 e 2 ( A -«)) 
+ (l~y 2 ) (G-l)e 2 ( A - 2 ") 
(l-y 2 )(G-l) 



(l-yl) (e 2 ^-l)+2y 2 ( e 2 ( A ^)-l) 



{ ' 2 ' (e 2(A " u) -iy 



^ 2 (<5? + Ql) + 2y Jl - y 2 Q Ze A - 2u (Q 2 - Q%) 



vl (Qj + Ql 
(l-y 2 e 2 ( A -^Y 



(6.30) 
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In order to simplify this formula we write, using (|6 . 14[) . (|6.16|) : 



(i - vl) 



2(A-2ti) 



- 1 



2yo (< 



2(A-u) 



-1 + (G-l)e 



2(A-2u) 



(1 - y 2 ) (1 - e~ 4 «) + 2y 2 (e 2 ^ - l) 
G 

[G + y 2 (l-G)] 
G 



= -(l-% 2 ) ^1 
+ (1 - y 2 ) (2Z 2 + Z 4 



[G + y 2 (l-G)] 



2y 2 



[G + y 2 (l-G)] 



- 1 



[G+ % 2 (i-G)]y 



( e a(A-«) _ ^ 



Z 2 + (l-G) (1 - yg) 



[G + y 2 (l-G)] ■ 
Plugging these formulas into (|6.30[) we obtain, after some computations: 



S _ yl{l-vl) (Qj + Qi) 
z 



(1 - y 2 e2(A-«)) 
Z 



(1 - y 2 ) (2Z + 



G 



[G + y 2 (l-G)] 



[G + y 2 (l-G)] 
l-vl)(G-l 



(l - y 2 e 2 ( A -«))' 

y 6 (Ql + Ql) 



Z {Q\ + Ql) + 2y a ^/l - y%e A - 2u (Q 2 - Q 
2(l-y 2 )Z(l-G) + Z 3 



(l - yge 2 ( A -")) [_ [G + y 2 (l-G)r 



(6.31) 



Summarizing, we have transformed the original problem (|3.15p . p. 161) . (|5 .4[) . 
(|5.11l) into the system of equations (16.24l) - (|6.27l) with A as in (|6.21l) . ^ as in 
(|OTj) . Ci as in and A, it given by ([H7T1J) . ([635)1 , The initial data for 

(Qi,Q 2 ,G,Z) are as in (jOSl) . 

Some of the forms that we have derived for the ODE problems above are more 
convenient for describing the solutions in different regions of the phase space. 
We will change freely between the different groups of equivalent variables in the 
following. 

6.2 Local existence of the curves 71, 72. 

With the reformulation of the problem obtained in the previous subsection the 
existence of the curves 71, 72 in a neighbourhood of the point (yo,Vo) can be 
obtained using standard ODE theory. 
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Proposition 7 For any yo £ (0, 1) and any j3o > there exist 6 > and two 
curves j%, 72 that can be parametrized as 



j i = {(y,V):y <y<y +S,V = V i (y)} , i = l,2 (6.32) 



with the functions V\ (y) , V2 (y) as m ( (5.5)) . d<5. 6]) satisfying \5.2\) , \5.3\) , 
the functions U, A satisfying \3.15)) , LS.16}) and \5.9j) . \5.1(ft) with p, p as in 
{5.1$ , i5.15\) and ax, a 2 solving i5.11\) . i5.1ty) . 

Proof. The arguments in Subsection 16.11 show that the proposition follows 
from proving local existence and uniqueness for (|6.24[) - (|6.27[) with initial data 
(|6.28[) . Since the right-hand side of (|6.24l) - (|6.27[) is analytic in a neighbourhood 
of (Qx, Q2, G, Z) = (1,1,1,0) it follows that there exists a unique solution of 
(H2H1), ([^Ml) - ([rT27)) on an interval of the form < x < S for some S > 0. 
Moreover, for such a solution A — > as \ — > + , whence Z ~ 2\ as \ — > + . 
Therefore (gl5| yields: 

s-2 X 2 as x~>0+ , as s^0+, 

Z ~ V2s as S4 + . (6.33) 

Using (|6.1[) it follows that: 

y-yo 



yo 



as y^yt (6.34) 



Combining then (jlTTj) and we obtain ([5T5]> . The asymptotics (pTTUf 
follows from the asymptotics for G, Z in an analogous way. ■ 

Moreover, we can prove Proposition [5] in a similar way. 

Proof of Proposition [11 It follows from (f3T25|) . (|6TTj) . (|Q3]) . (J04J. ■ 

We notice for further reference that we have also proved the following result: 



Proposition 8 There exists a unique solution of the system h6.S\) - ^75]) with Q 
as in ifi.ty) and initial data {Qx, Q2, A, u) = (1, 1, 0, 0) as s — > + . 

6.3 Steady states for the system ( 1Q41) - ( 16^271 ). 

In order to study the steady states of (|6.24l) - (|6.27[) it is more convenient to use 
the form of the equations in (|6.2[) - (|6.5[) . Then the steady states are characterized 
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by: 



QiCi = i = 1,2, 



(6.35) 



-2A 



1 = - 

2 



+ 



Qi 



Cie" + y e A VCf+l 



[d 2 + 1] 



Ql 



«2 



3e 



-2A 



1= 2 



8 (CO 2 



?I (C 2 ) 



Cie u + j/oeVCi +1 Cae" + y e A \/C|Tl 
The first and third equations imply: 

3e~ 2A -1 = 0. 
Then, the second equation reduces to: 

Qi Q\ 



(6.36) 
(6.37) 

(6.38) 
(6.39) 



Notice that (|6.39p implies that at least one of the variables Qi, Q2 is different 
from zero at the steady state. Suppose that both of them are different from 
zero. Then £1 = £2 = 0, whence, using 



Cf + 1 + J/oCie 



1 



Vl 



1,2 , Ci <C 2 



it follows that: 

and (I6.39[) reduces to: 



?/5 



(6.40) 



(Qi + Qi) = 



Ay e A iy V3 



39 39 

This defines a family of steady states. Local analysis near these solutions 
indicates that they are reached for finite values of y. Since we are interested in 
solutions defined for arbitrarily large values of y > yo a more detailed analysis 
of these solutions will not be pursued here. We will then restrict our analysis 
to the solutions for which Q1Q2 = 0. 

Suppose that Qi ^ 0. Then d = 0. (f6~4"0"]) implies: 



C2 



1. 



Voe J 



1 - \/C 



< 0. 
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This contradicts d < Therefore for solutions with Q1Q2 = we must have 
Q 2 ^ whence C2 = 0. Then (jOOf is satisfied and (|6~3T)| yields: 



/4V3 



12 = 



■yo 



2^ 



We remark that for this solution: 



2/ie A ° 



2y n V3(l-y 2 ) 



In order to have Ci < C2 = we need yo £ (0, |) . 

Summarizing, for each yo £ (0, 5) the system (|6.24p - (l6.27[) has the following 
steady state: 



Q2 = Q2.00 : 
A = Aoo = 



= 0, 



2^ 



log (3) 



u = U Q 



log I y 1 — i/g 

We also introduce the following notation for further reference: 

_ 2fre A ° 

0. 



(6.41) 
(6.42) 

(6.43) 
(6.44) 



Cl.oo 
C2,oo 



(/l2_e2Ao.) 



2yoV3(l-yg) 
l-4y 2 ' 



(6.45) 
(6.46) 



6.4 Linearization near the equilibrium. 

The main result that we prove in this subsection is the following: 

Theorem 9 For each yo £ (0, V\ the point Pi = (Qi >0 o, Q2,oo, Aoo, u<x,) defined 
by ft6.41\ )- f6-44\i * s arl unstable hyperbolic point of the system H6.2\) - [U75\) . The 
corresponding stable manifold of the point (Qi }0 o, Q2,oo, Aqq, Uoo) that will be 
denoted by A4e is three-dimensional and it is tangent at this point to the subspace 
generated by the vectors 

















V J 





( 





'— — 3 — 

~ 3 







■5 — 3 — 

1 



> . 



/J 



(6.47) 
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Proof. The key ingredient in the proof of this theorem is the linearization 
of the system (|6.2p - (|6.5|) around the point (Qi t oo, Q2,oo Aqq, Uqo) . Let us write: 

A = Aoc + L, 
u = u oa + v, 

Qi = Qi,oo + qi = qi, 

Q 2 = Q 2 ,oo + 92- 

Neglecting terms quadratic in \L\ + \v\ + \qi\ + \q 2 \ we obtain, after some tedious, 
but mechanical computations, the following linearized problem: 



dgi 
ds 

dq 2 
ds 



2h 2 



2 1 



(V - 3) 
2 (1 ~ Vl) 

■qi -+ 



(1-4^) 



Vl) 

-57«l, 



;i - vl) 

3W§ 



= 3£. 



(6.48) 
(6.49) 

(6.50) 
(6.51) 



Looking for solutions of the linearized problem with the form: 

( A x \ 

A 2 
A 3 
V A, ) 

we obtain the following possible values of 7 with their corresponding eigenvec- 
tors: 



2 1 



7i 



y 2 ) 



(A 1 \ 

A 2 
A 3 
V ^ J 



( 1 \ 





V J 



72 



-2 o 



/ Al \ 




A 2 




A 3 









\ 

(1-fo) 



73 



MO 








A 3 













5 3 

_VEA 

1 



33 



74 



yo 





/ 





\ 














5 3 






















V 


1 


/ 



The theorem then follows from standard results for stable manifolds (cf. for 
instance [4], [19]). ■ 

6.5 Reformulation of the solution in the original variables. 



Our goal now is to obtain a trajectory connecting the point (Qi, Qi, A, u) = 
(1,1,0,0) at s = with the point Pi at s — oo for a suitable value of 9 (or 
equivalently (3o). Let us remark that such a trajectory would satisfy the re- 
quirements in Theorem Q] Indeed, notice that such a trajectory behaves near 
the point (yo,Vo) as stated in Theorem Q] due to Proposition [71 On the other 
hand, such a trajectory would belong to the stable manifold of the point P\ and 
therefore its asymptotic behaviour as s — > oo would be given by: 



( Qi \ ( Qi.oo \ 



A 

u 



V u J 



\ 



Aoo 



-2s 



( o \ 



2 
3 
-1 



2 ( 1 -yp) 
C 2 e t 1 " 4 ^) ' 



/ 1 \ 





V o J 



for sufficiently small yo (cf- [I])- Notice that the smallness of yo guarantees that 
the last term yields a contribution larger for ,s — > oo than the first quadratic 
corrections if C-2 ^ 0. 

Using (|6.1j) we obtain the following asymptotics for the original set of vari- 
ables U, A, (Ti , Vi, i = 1, 2 : 



A log ^V3^ as 



y -)• oo, 



U 2 



C 2 



?7o 



as y — > oo, 



Q2,oo — as y^oo, 



2W3(l-yg) 

(i-4% 2 )y 



+ o (1) as j/4 oo, 



as y — > oo, 



V 71 -2/o f ^0 



V3yo y \y 



— I as y — > oo. 
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in particular these formulas prove Theorem [5] 



6.6 The shooting argument: Approximation of the stable 
manifold Aig for small yo. 



Since the stable manifold Me is three-dimensional we cannot expect the point 
{QiiQi, A, u) = (1, 1, 0, 0) to belong to Me for generic values of 9. The intuitive 
idea of the proof which follows is to show that the manifold Me divides the set 
{0 < G < 1, Z > 0, Qi > , i = l,2} into two different regions. If the point 
(1,1,0,0) lies on different sides of Me for different values of 9 then by continuity 
there must exist a value 9* of 9 such that (1, 1,0,0) G Me- In the rest of the 
paper we will obtain approximations to the manifold Me for yo small that will 
show that the point (1, 1,0, 0) lies on different sides of Me for large positive 
values of 9 and small positive values of 9. More precisely, the main result of this 
subsection is the following: 

Theorem 10 There exists yo small enough such that, for any yo in the interval 
[0,j/o] there exists 9* = 9* (y ) > such that (1, 1,0,0) G Me*- 

Proof. In order to prove Theorem [TU] it is convenient to use the coordinates 
(Qi,Q2,G, Z) (cf. (|6.9[) . (I6.14P 1 ). These variables satisfy the system of equations 
(|6.24[) - (|6.2T[) . The steady state P\ — P\ (yo) is given in these coordinates by: 



■Pi = (<9l,oo) Q2,oa, Goo, Zoo) = ( 0, —3-^ 



3iV£'3'V (1-l/g) 



(6.52) 



The point Pi depends continuously on yo if yo G [0, |] . If yo = the system 
(j6TM]) - (|6T2Tj) becomes: 



dQi 
dQ 2 



dG 



dZ 



-2GZQ U 
-2GZQ 2 , 



—— = 2G 



Z(l-G) 



3G- 1 



(Qi + Qi) 



[z^ + iy 



■z (Qi + Qi) )(z 2 + i 



(6.53) 
(6.54) 

(6.55) 

(6.56) 



Theorem [9] shows that the point Pi (yo) is hyperbolic for yo G (0, |] with a 
three-dimensional stable manifold Me = Me (yo) ■ On the other hand two of 
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the eigenvalues associated to the linearization around Pi of the system (|6.24l) - 
(|6.27l) degenerate for y — 0. More precisely, let us write G = \ + g- Since 
Pi (0) = (0,0, i,0) we obtain the following linearization of (|6.53[) - (|6.56[) near 
Pi(0): 

dQi _ n dQ 2 _ dG _4Z dZ 
dC ' d( ' d( ~ 9 ' dC " 5 ' 



The corresponding eigenvalues are |o, 0,— ^fii^fi^ an d the corresponding 

/1\ /o\ / o \ f o \ 
1 o 

1\[Z 



eigenvectors are 





V o J 







v i y 








Standard results 



v i y 



(cf. [4]) show the existence of a centre-stable manifold that will be denoted by 
A4e (0) that is invariant under the flow defined by the system f|6.53[) - (|6.56[) and 

f/ X \ f°\ ( ° \ 

is tangent at Pi (0) to the plane spanned by 






v»y 



i 



v o y 





2x/3 



v i y 



Classical results (cf. [4]) then show that it is possible to obtain a con- 



tinuously diffcrcntiable four-dimensional manifold Mg,cxt 
(Uo, Qi,Q2, G, Z) e M ex t such that: 

M e ,ext n {y = b} = M 9 (6) 



c [o,|] 



with 



(6.57) 



for any b G (0, |) . Indeed, the manifold A4e, cx t is any centre-stable manifold 
at the point (yo, Qi, Qi, G, Z) = (0,0, 0,|,0) associated to the system (|6.24l) - 
(|6.27l) complemented with the additional equation 



dyo 
d( 



= 0. 



(6.58) 



More precisely, we make use of the fact that the dynamical system of interest 
has a smooth extension to an open neighbourhood of the stationary point un- 
der consideration. The manifold A^e,ext is the intersection of a centre-stable 
manifold for the extended system with the subset defined by the inequality 
J/o ^ 0. The manifold A4 cx t contains all the points of the form (yo,Pi (yo)) 
with yo € [0, 5] since they remain in a neighbourhood of (0, 0,0, |,0J for ar- 
bitrary times. Moreover, the manifolds Ate, ex t H {yo = b} are invariant under 
the flow (|6.24[) - (|6.27[) and since they are formed by points that remain in a 
neighbourhood of (0,0,0, |,0) for arbitrarily long times, it follows from (|6.58[) 
that the points in Me.cxt H {ya — b} are contained in the stable manifold asso- 
ciated to the point Pi (yo) ■ The uniqueness of the stable manifold then implies 
Me (b) C Me,cxt n {yo — b} . Moreover, the form of the tangent space to Mg, e xt 
at the point (0,0, 0, |,0) implies that the dimension of Me, e xt H {yo = b} is 
three for small b. Since this is also the dimension of ftAg (b) the relation (|6.57p 
follows. The continuity of -M C xt then implies that the centre-stable manifold 
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Mg (0) can be uniquely obtained as limit of the manifolds Mg (yo) as yo — ^ + . 
In particular the manifold Mg (0) is unique. 

The properties of the manifold Mg (0) can be analysed in more detail. We 
remark that the curve: 

^(Z 2 + 1)VG(1-G) = 4- . Qi = Q 2 = (6.59) 

belongs to Mg (0) since the hyperplane {Qi = Q2 = 0} is invariant under the 
dynamics induced by (|6.53p - (|6.56p . On the other hand, the invariance of (|6.53p ~ 
(|6 . 56|) under rotations in the (Qi, Q2) -plane allows the problem to be reduced 

2\ 



to one with smaller dimensionality. More precisely, defining <3 = y 5 [Q\ + Q 
leads to the system: 

^ = -2GZQ, (6.60) 
dC r 3 1 

-^ = 2G Z{1-G)-6[Z 2 + 1Y Q 2 , (6.61) 

^ = (3G - 1 - 9ZQ 2 ^{Z* + 1)) (Z 2 + 1) . (6.62) 

We will denote by Afg the (two-dimensional) invariant manifold associated to 
the system (|6.60p - (|6.62p that is obtained from Mg by taking the quotient by 
rotations in the Qi and which contains the curve (|6.59p . 

Our goal is to show the existence for any y sufficiently small of a value 
9* = 9* (yo) of 9 such that the manifold Mg* (yo) contains the point Qi — 
Q2 = 1, G = 1, Z = 0. This will be done by showing that the corresponding 
statement holds in the case yo = and then doing a perturbation argument. 
The statement about the manifold A4g* (0) is equivalent to the statement that 
Afg* contains the point (1, 1,0). It will be shown that the latter statement is 
true and, moreover, that when 9 is varied through the value 9* the manifold 
Afg moves through (1, 1,0) with non-zero velocity. It then follows that 7^(0) 
moves through (1, 1,1,0) with non-zero velocity. Note that the coefficients of the 
system extend smoothly to an open neighbourhood of the manifold Mg* (0). As 
a consequence the manifold A4g, cx t extends smoothly to small negative values of 
yo- The desired statement concerning A^e(yo) is a consequence of the implicit 
function theorem. In more detail, the statement that Mg depends on 9 and yo 
in a way which is continuously differentiable means that there is a C 1 mapping 

from the product of a neighbourhood of (0,6**) in W 2 with Mg(0) into a 
neighbourhood of (1, 1, 1,0) with the properties that its restriction to yo = 
and 9 — 9* is the identity and that the image of {(yo,@)} x Mg*(0) under 
^ is Mg(yo). The condition that the manifold moves with non-zero velocity 
implies that if xq denotes the point of Mg*(0) with coordinates (1, 1, 1,0) the 
linearization of \1/ at the point (0, 9* , xo) with respect to the last four variables 
is an isomorphism. This allows the implicit function theorem to be applied. 

In order to check the existence of 9* it is enough to study the behaviour of 
the manifolds Afg for 9 — > + and 9 — > 00. These manifolds are two-dimensional 
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manifolds in the three-dimensional space (Q, G, Z) . Notice that the structure 
of the manifolds Afe can be easily understood using the fact that the parameter 
9 can be rescaled out of the system (|6.60[) - (|6.62p using the change of variables: 



Then (|QD|l - ([6T62]l becomes: 

-2GZq, 



dq 

d( 
dG 

dz 



2G 



Z(l-G)- [Z 2 + 1]' 



3G - 1 - + [Z 2 + 1) 



(6.63) 

(6.64) 
(6.65) 
(6.66) 



Let us denote by Af the centre-stable manifold at the point (q, G, Z) = 
(0, §,0) for the dyn amies (I6.64[) - (I6.66[) . The manifold Af contains the curve 

{(Z 2 + l)G(l-G) 2 -^, q = o\- Notice that: 



V8Q,G,Z) e Af 



(Q,G,Z)eAf e 



Therefore the family of manifolds Afe can be obtained from the manifold Af 
by means of the rescaling (|6.63|) while keeping the same value of the variables 
G, Z. In order to check if (Q,G,Z) = (1,1,0) 6 Afe we just need to describe 
in detail the intersection of the manifold Af with the line {G = 1, Z = 0} . Once 
the existence of a value 9* of 9 for which the manifold Afe* contains the point 
(1, 1, 0) has been shown the statement that the manifold Afe moves through this 
point with non-zero velocity follows immediately from the rescaling property. 

Notice that the plane {q = 0} is invariant for the system of equations (|6.64p ~ 
(|6.66p . The analysis of the trajectories of (I6.64|) - (I6.66|) in this plane can be done 
using phase portrait arguments. There is a unique equilibrium point at (G, Z) — 

(§, 0) with stable manifold j (Z 2 + l) G (1 - G) 2 = ^-j . This manifold splits 

the plane (G, Z) in two connected regions. The trajectories starting their motion 
in the region that contains the point (G, Z) — (0, 0) reach the line Z = for 
a finite value of £ if Z > initially and eventually develop a singularity where 
Z approaches — oo at a finite value of £. On the other hand, the trajectories 
starting their motion in the region containing the point (G, Z) = (1, 0) move in 
the direction of increasing Z towards Z = oo, G = 1, a value that is achieved 
for a finite value of £• 

Notice that the solutions of (I6.64[) - (I6.66[) starting their dynamics in the set 
{0 < G < 1, Z > 0} can only evolve in two different ways. Either the trajectory 
remains in the region where Z > for arbitrarily large values of £ or the 
trajectory enters the region {Z < 0} . In the second case this can only happen 
through the set G < |. Since G is decreasing it remains in the set {Z < 0} for 
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larger values of £ and eventually it approaches Z = — oo for some finite value of 

c 

Suppose otherwise that the trajectory remains in the region where Z > for 
arbitrary values of Q. Then q decreases to zero and the behaviour of the trajec- 
tories is then similar to the ones in the plane {q — 0} . We now claim that either 
this trajectory belongs to the stable manifold JV or it satisfies limf_>.f » Z (£) = oo 
for some < oo. In order to avoid breaking the continuity of the argument we 
will prove this result in Lemma [TT] in Section [7] 

We will show that there exists a point of the line {G = 1, Z = 0} in the 
manifold JV. The points of this line enter the region {0 < G < 1, Z > 0} due 
to the form of the vector field associated to (|6.64[) - (|6.66p . If q (0) > is small, 
Lemma [T^l shows that Z approaches Z — oo for a finite value of £. Suppose now 
that q (0) > is sufficiently large. Then the trajectory enters the region {Z < 0} 
for a finite value of ( as the following argument shows. A solution which starts 
at (§0)1)0) with go large immediately enters the region Z > 0, G < 1. The 
inequality Z < 1 will hold for at least a time i since 4? < 4 as long as Z < 1. 
The aim is to show that for g sufficiently large Z will become negative within 
the interval [0,4]. From now on only that interval is considered. Integrating 
the equation for q gives the inequality g(£) > e~?qo. The equation for G then 
shows that G(C) < e -0 '* 1 ^ where a(qo) = q^e^ 1 — 1. Choose go large enough so 
that e -TO a (<ro) < i. when C = jo the inequality Z < ± still holds. Under the 
given circumstances G is decreasing on the whole interval [0, j}. The equation 
for Z shows that by the time C = ^ at the latest Z has reached zero. 

Let Ui be the set of positive real numbers go for which the solution starting 
at (go, 1,0) is such that Z — > — oo as ( — > (* , where denotes the maximal 
time of existence, and let U2 be the set of positive real numbers go for which 
the solution starting at (go, 1,0) is such that Z — > +00 as £ — > £*. It follows 
from Lemma fl~3l that U2 is open. We also know that U\ is open. Moreover, it 
has been proved that both U\ and U2 are non-empty. By connectedness of the 
interval (0, 00) it follows that there must be a value of go for which the solution 
starting at (go, 1,0) is neither in U\ or Ui- For that solution Z is non-negative 
and does not tend to infinity and thus, by Lemma [Til it is the desired solution 
which lies on JV. 

The equivalence between the existence of the self-similar solution described 
in Section[5]and the existence of a trajectory connecting the points (Qi, Q2, G, Z) = 
(1,1,1,0) and (Q1.00, Q2,oo) Aqo, Uoo) proved in Subsection 16.51 concludes the 
proof of Theorem [TU1 Theorem [JJ is just a Corollary of Theorem [TD1 ■ 



7 Some auxiliary lemmas used in the analysis of 

(ESlD-dEBSD. 

Lemma 11 Suppose that a solution of \6.6Ji\ )- $6. 66\) is defined for £ £ [C*> C*) ) 
where is the maximal time of existence. Suppose that Z (£) > for any 
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C G [£*,C*) an d a l so that G(C*) G (0,1), <z(C*) > 0- Then, either the curve 
{(q (£) , G (£) , Z (0) :£G (C*,C*)} * s contained in the stable manifold J\f or 
lim£_>.£* Z (£) = oo. 

Proof. The plane {G — 0} is invariant under the flow associated to (|6.64l) - 
(|6.66l) . On the other hand, the vector field on the right-hand side of (|6.64p - (|6.66l) 
points into the region {G < 1} if q ^ 0. Therefore the region {0 < G < 1, q > 0} 
is invariant for the flow defined by (|6.64[) - (|6.66p and we can assume that the 
inequalities < G (£) < 1, q (Q > hold for any ( e [C*> C*) • We now have two 
possibilities: 

limsupZ(C) < oo , (7.1) 
limsupZ(C) = oo. (7.2) 

Suppose first that (|7.ip holds. Then, there exists M > such that 

Z(Q<M for any Ce[C.O- ( 7 - 3 ) 

We claim that in this case the trajectory {(q (£) , G (C) , Z (£)) : £ g (£*,£*)} i s 
contained in TV". Notice that in this case, the boundedness of \(q, G, Z)| implies 
that C* = oo. Since {GZq) (C) > for ( € [C*,oo) it follows from ([fTM]) that 
q (C) is decreasing. Therefore qoo = lim^oo q (£) exists and is non-negative. 
Suppose that < goo. Then < qao < q (() for any ( e [C*,°°) • Integrating 
(|6.64p we obtain (GZq) (() d( < oo, whence 

(GZ)(C)dC<oo. (7.4) 

Since ^, ^ are bounded, CO} implies lim^oo (GZ) (£) = 0. Then (j^5l) 
implies: 

^<VG 
rfC " 00 

for C > Co sufficiently large. Therefore lim ( ;_ i . 00 G (£) = 0. Equation (|6 . 66[) then 
yields: 

dZ 1 
~d( - ~2 

for £ > Co large enough. Then Z (£) < for large £, but this contradicts the 
hypothesis of the lemma. It then follows that qoo = 0. 

Due to (|7.3p and since lim^oo q (£) = we can approximate the trajec- 
tories associated to (|6.64|) - (|6.66[) for large values of C using the corresponding 
trajectories associated to fl6.64[) - (I6.66|) for q = 0. The study of the trajecto- 
ries associated to (|6.64p - (|6.66p that are contained in {q — 0} n {0 < G < 1} 
reduces to a two-dimensional phase portrait. These trajectories can have only 
three different behaviours. Either they are contained in Af H {q = 0} , or they 
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reach the plane {Z = 0} , with G < |, entering {Z < 0} , or they become un- 
bounded. The continuous dependence of the trajectories with respect to the 
initial values as well as the fact that lim^oo q (£) = implies then that either 

lim^oodist ((?(C),G (£),Z(C) ),Mn{q = 0}) = 0, or 2 (C) < for some 
£ < oo, or Z (() > M + 1 for some £ < oo. The second alternative con- 
tradicts the hypothesis of the lemma. The third alternative contradicts (|7.3|) 
and therefore only the first alternative is left. However, in that case lim.f_>oo 
(q (C) , G (£) , Z (C)) = (0, ^, 0) and the trajectory is contained in Af as claimed. 

Suppose then that (I7.2[) holds. We claim that in this case lim£_>.£* Z (£) = oo. 
Notice that the monotonicity of q (£) implies that limj_j.^* q (£) = q^ exists. 
We will first prove that q^ — 0. Suppose that, on the contrary, q^ > 0. Then 
9 (C) > loo > for any ( £ [C* , C*) ■ Equation (|6.66[) as well as G < 1 yields: 

dZ 



( , c < (2-ZqlVZ^+T) (Z 2 + l) 

for any ( e [C*>C*)- This inequality implies ^5 < for Z > Z^ = Z^ (g^) . 
Therefore Z (£) < for ( 6 [C*>C*) an d this contradicts (|7.2I) . From now on 
take goo = 0. We can then assume (|T.2[) and 



hm g (C) = 0. (7.5) 

Suppose also that liminf^-^. Z (£) < 00. This is equivalent to the existence of 
< M < 00 and a subsequence {£„} with £„ — > (* as n — > 00 such that: 

Z(C«)<M. (7.6) 

We now claim that: 

UmjZ(C)?(0]=0. (7.7) 
To prove (|7.7p we argue as follows. Combining (|6.64[) , (|6.66l) we obtain: 
d 



(Zq) = qZ 2 (G-l) + q (3G - 1) - Zq 3 yj (Z* + if. (7.8) 



We now use the inequality Zy {Z 2 + 1) > Z 4 for Z > 0. Then, using also the 
inequality G < 1 : 

^-(Zq)<q- 1 \(3G-l)q 2 -(Zqf}. (7.9) 

It follows from this inequality, as well as (I7.5[) that for any e > 0, every tra- 
jectory satisfying the hypothesis of Lemma [TT] and entering any of the regions 
{(q, G, Z) : Zq < e } for £ sufficiently close to £* remains in such a region for 
later times. If (* = 00, the meaning of sufficiently close is large enough. Due 
to (|7.5p and (|7.6p . for any e > 0, there exist ( n arbitrarily close to (* such that 
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(Zq) (£„) < e. Then (Zq) (Q < e for any ( 6 (CmC*) • Since e is arbitrary we 
obtain (|7Jl) . 

Combining (|7.5p and (|7.7p it follows that: 

lim^i (C) = Km 6 2 (C) = (7.10) 

where: 

We can then rewrite (|6.64p . (16.66)) as: 

§ = -2GZg, (7.11) 

— = 2G[Z(1-G)-(Z + 1)6 2 (0], (7.12) 

^ = (3G-1-<5 1 (C))(Z 2 + 1). (7.13) 

We now claim the following. Given any £0 belonging to the interval (0, |) 
suppose that the trajectory under consideration enters the set: 

£l ea = Ig> i+e , Z> 1 

for some C < C* sufficiently close to (* . Then lirrif_>f* Z (() = 00 and < 00. 
The proof as the follows. Due to (|7.10[) the set fl £o is invariant for (|7.11[) - (|7.1 3[) 
if C is close to (*. Then, for £ close to (* we have: 

f + 

and this implies lim^_^^» Z (£) =00 and £* < 00. 

Therefore, to complete the proof of Lemma [TT] it only remains to prove that 
the trajectory enters fi eo for values of C sufficiently close to C* ■ Due to (|7.2p and 
(|7.6p there exists a sequence {Cn} with ( n < (n < C* such that: 

Z (C„) = 2M and — (Q > 0. 



Due to (|7.13p this implies: 



lim sup G(C„) > \. (7.14) 

n— >oo o 



On the other hand, a Gronwall type of argument applied to (|7.13p implies the 
existence of olm > 0, depending only on M such that: 

< y < Z (C) < 4M for C S [Cn, Cn + a M ] . (7.15) 
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Comparing the solution of the equation 1 2[) with the solution of the equation 
^ = 2GZ (1 — G) with the same initial datum at £ = Q n and taking into ac- 
count (|7.14p . (|7.15l) it then follows that, for n large enough [q (( n ) , G (£„) , Z (Cn)) S 
f2 £o . Therefore lim£_^» Z (() = oo. This contradicts (|7.6I) and the lemma follows. 
■ 

Lemma 12 There exists 5 > sufficiently small such that, the solution of 
H6.64\ )- (6- 66]) with initial value (q (0) , G (0) , Z (0)) = (<jo, L 0) and < qo < 
satisfies lim^^. Z (£) = oo, where C* denotes the maximal time of existence of 
the trajectory. 

Proof. The trajectory enters the region {Z > 0} and as long as it remains 
there, the function q (£) is decreasing. The inequality ^ < 4 holds as long 

as Z < 1. It follows that Z < 1 on the interval [0,^1. On that interval the 
inequality dflgsG) > — 2^q^ holds and hence G > e~ q ° . Furthermore 

|l>3e-«8 -l-V2&=P(q a ). (7.16) 

Choose S sufficiently small that (3(6) > 1 and e~ s > |. Then Z (i) > i and 
G > i on [0, |] . Choose e > and suppose that 2<5 2 < e 4 . Then it follows from 
(|7.9I) that the set defined by the inequality Zq < e is invariant. Thus the solution 
remains in that region on its whole interval of existence. Now 5i(C) < eV e 2 + S 2 
and ^(C) < (e 2 + 5 2 ). Let [0, Ci) be the longest interval on which G > 4. From 
what has been shown already (i > \, Reduce the size of e if necessary so that 
eVe 2 + S 2 < \. Then it follows from (|7.13[) that Z is increasing on [0, Ci) and 
hence is greater than \ for ( > Ci- Putting this information into (17.12)) shows 
that provided e 2 + q 2 < then G cannot decrease. For 5 sufficiently small 
this gives a contradiction unless £i = C* ■ I R particular there is a positive lower 
bound for Z at late times. Furthermore 1 3[) implies that lim^_y^* Z (£) = oo 
and the lemma follows. ■ 

Lemma 13 Suppose that a solution satisfying the hypotheses of Lemma \ll\ with 
C* = has the property that lim^_j.^* Z(Q = oo. Then any solution starting 
sufficiently close to the given solution for £ — also has the property that Z 
tends to infinity on its maximal interval of existence. 

Proof. To start with a number of further consequences of the hypotheses of 
Lemma [TT1 will be derived. The assumption on the initial condition only plays 
a role towards the end of the proof. It has been shown in the proof of Lemma 
[TT1 that lim^-^. q (Q — 0. We now claim that (|7.7p holds. Suppose that it is 
not true. Then we claim that the limit lim^-^* (Zq) (£) = L exists and that 
L > 0. Indeed, notice first that liminf^-j.^* (Zq) (£) > 0. Otherwise there would 
exist a sequence {Cn} such that lim n _ i . 00 £„ = (* with limn^oo (Zq) (£„) = 0. 
Combining this with the fact that q — > and (|T.9[) we would obtain (|7.7p . a 
contradiction. Thus liminf^^* (Zq) (C) > 0. Using again the fact that q — >• 
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and (|7.9[) it follows that (Zq) is monotone decreasing for £ close to C*, whence 
the limit liirif_>.f» (Zq) (£) = L exists. Moreover we have obtained also in this 
case that (Zq) (£) > L for ( close to (* ■ 

It follows from the proof of Lemma HT1 that £* < oo. By the boundedness of 
the right hand side of (|6.64p it follows by integrating this equation between £ and 
(* that q(C) < a^ 1 ((*-() for a positive constant a. Hence q _1 (C) > a(C* — C) -1 - 
This can be used together with the limiting behaviour of Zq to estimate the 
right hand side of (j7.8[) from above. The first term is negative and can be 
discarded. The second term tends to zero as ( — > £*. The third term can be 
written in a suggestive form as — q~ 1 [(Zq) ^J((Zq) 2 + q 2 ) 3 ]. The expression in 
square brackets tends to a positive limit as £ — > (* . Thus the right hand side 
of (|7.8|) fails to be integrable, contradicting the fact that Zq is positive. This 
contradiction completes the proof that lim£_>£* (Zq) (() = 0. 

We now use some arguments analogous to the ones used in the proof of 
Lemmas [TT] and [HJ As a next step we prove that G(£) tends to a limit as 
£ — > £* and that this limit is greater than |. We first claim that: 

S = limsupG(C) > -. (7.17) 

Indeed, suppose first that S — limsup^^. G (() < |. Since lirri£_>.£. (Zq) (Q = 
we can approximate ([51)D]) - ([51j2]) by the system (fTTTTT) - (f7TT5]) . Using (fTTT^l) it 
follows that Z (() is decreasing for £ close to (* . This contradicts (|7.2[) and then 
(|7.17p follows. On the other hand (|7.12[) implies that G is increasing if G > \ 
for £ close to (* . Using (|T. 1 7[) it then follows that G increases for £ close to £* . 
Therefore the limit lim^_j.^* G (£) exists and: 

c lim G(C )4 

Since G is monotonically increasing we can parametrize Z as a function of 
G. Let us denote the corresponding function by Z = Z (G) . Then by 1 2[) and 
CLHD: 

jOggj) = (3G-i-^(Q)(i + ^- 2 ) . lg) 

dG 2G[(1 - G) - (1 + Z-^)S 2 (C)}' 

If the limit of G were less than one the right hand side of this expression would 
be bounded and it would follow that Z was bounded, a contradiction. Hence 
lim c ^.G(C) = l. 

To complete the proof the condition on the initial data in the hypothe- 
ses of the lemma will be used. Since lim^-^* Z (£) = oo, lim£_>£* q (£) = 0, 
lim^-^* (Zq) (£) = and lim^-^* G (Q > -| it follows that for any sufficiently 
small S > and for any solution (q, G, Z) that is sufficiently close to (q, G, Z) 
at £ = we have for some Co < C* : 

q((o)<5 3 , G(Co)>^+<5 , (Zg)(Co)<* , ^ (Co) > |- (7.19) 
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It will now be shown that for 6 sufficiently small the region defined by these 
four inequalities is invariant. On the part of the boundary of the region where 
q = 5 3 we have || < 0. On the part of the boundary where Zq — S assuming 

that 5 < suffices to show, using (|7.9[) . that the derivative of Zq is negative. 
On the part with G = \ + S the following inequality holds: 



dG 2 
~d( ~ 3 



2 

35- 1 



(Si + 5*)* - (S + S 3 ) 2 



(7.20) 



Choosing 5 sufficiently small implies that the right hand side of this inequality 
is positive. On the whole region 



V / S 2 + S 6 ). 



(7.21) 



If 6 is small enough then this quantity is positive. Putting these facts together 
shows that the solution starts in the region of interest when C = Co and stays 
there. In particular G (() > -| + 5 for ( > (o- Therefore Z blows up in finite 
time due to (|7.13j) and Lemma IT3l follows. ■ 
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